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1. Introduction 

The synthesis of superheavy chemical elementJi"^] jn the laboratory is one of the 
biggest challenges in nuclear physics. It is an attempt for clarifying the existence 
limits of all the chemical elements, as well as the completion attempt of the periodic 
table of chemical elements. We are concerned with heavy-ion collisions 

248Cm + 48Ca 

with different impact parameters in this paper. Let A and Z be the mass number and 
the proton number, respectively. The neutron number N is defined hy N ~ A ~ 
so that N/Z of ^^*Cm and **Ca are 1.58 and 1.40, respectively. If fusion appears, 
296lv 296^15) with N/Z = 1.55 is produced. 

Fast charge equilibratiorP is expected to appear in low-energy heavy-ion reac- 
tions with an incident energy of a few MeV per nucleon. It provides a very strict 
limitation for the synthesis of superheavy elements. Actually, the N/Z of the final 
product is not above nor below the N/Z of the projectile and the target (1.40 < 
N/Z < 1.58 in this case) in the case of charge equilibration, so that the proton- 
richness of the final product follows. Although the actual value of N/Z depends on 
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the two colliding ions, its value for the merged nucleus tends to be rather proton- 
rich for a given proton number of the merged system. This feature is qualitatively 
understood by the discrepancy between the /3-stability line and the N — Z-line for 
heavier cases. Superheavy compound nuclei are very fragile and fission is a very 
frequent channel which leads to disintegration of the compound nuclei even at low 
excitation energies. In this paper, following the evaporation prescription shown in 
Ref.,'^the fissibility (fission probability) of compound nuclei is calculated based on 
the time-dependent density functional theory (TDDFT). The presented method 
brings about the microscopic derivation of the fissibility. For reference, a possible 
fission dynamics based on the evaporation prescription was demonstrated in 

2. Methods 

2.1. Treatment of the thermal property 

Self-consistent time-dependent density functional calculations are employed in this 
paper. TDDFT reproduces the quantum transportation due to the collective dy- 
namics. In this sense, what is calculated by the TDDFT can be regarded as products 
after several 10~^^ s, which corresponds to a typical time-scale of low-energy heavy- 
ion reactions (1000 fm/c), as well as to the inclusive time interval of any collective 
oscillations such as giant dipole resonance, giant quadrupole resonance and so on. 
Meanwhile, thermal properties such as the thermal instability are not directly taken 
into account in TDDFT. Indeed, the Skyrme type interaction used in TDDFT (for 
example, see Ref.'^ is determined only from several densities. It is important that 
the most effective cooling effect arises from the emission of particles, and therefore 
it is expected that the break-up or fission of fragments including rather high inter- 
nal excitation energy is suppressed in the TDDFT final products. The additional 
thermal effects leading to the break-ups of fragments should be introduced. 

Here is a fact that simplifies the treatment of thermal effects, that is, the differ- 
ence of the time-scales. Different from the typical time scale of low-energy heavy-ion 
reactions, the typical time-scale of the thermal effects is estimated by the typical 
time interval of collision- fission (fission appearing in heavy- ion collisions): several 
10""'^^ s. It is reasonable to introduce an evaporation prescription simply to the 
TDDFT final products. In this context the TDDFT final fragments have the mean- 
ing of products just after the early stage of heavy-ion reactions (several 10~^^ s). 

2.2. Evaporation prescription 

In complete fusion reactions the cross-section for the formation of a certain evapo- 
ration residue is usually given by three factors!^ 

(^ERiEcm) = X! '^CPiEcm, J) X PcN{Ec,n, J) X PsviEcn, J) (1) 
J 

where ctcp, Pcn and Psv mean the capture cross-section, the probability for the 
compound nucleus formation, and the probability for survival of the compound 
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nucleus against fission. All three factors arc functions of the ccntrc-of-mass energy 
Ecm and the total angular momentum J, where J can be related with the impact 
parameter. For light systems Pcn and Psv are about unity and asR ~ J^j'^cp- 
But in superheavy systems the strong Coulomb repulsion and large angular mo- 
menta lead to small values of Pqn and Psv and therefore to the small cross-sections 
of the evaporation residues observed in the experiments. This means, different from 
light systems, it is necessary to introduce additional thermal effects for the super- 
heavy element synthesis. First, acp is sufficiently considered in the TDDFT if we 
restrict ourselves to a sufRciently high energy exceeding the Coulomb barrier (cf. 
sub-barrier effects such as tunnelling are not taken into account in the TDDFT). 
Second, Pcn is fully considered in the TDDFT, which is a kind of mass equilibra- 
tion also related to charge equilibration. Third, Psv whose relative time-scale is 
by no means equal to the former two probabilities is not satisfactorily considered 
in the TDDFT. This probability is much more related to thermal effects. Conse- 
quently, further consideration is necessary only for Psv as far as the energy above 
the Coulomb barrier is concerned. 

Several factors are included in Psv such as probabilities for fission of the com- 
pound nucleus, neutron-evaporation, proton-evaporation, deuteron-evaporation, 
alpha-particle-evaporation and so on. Probabilities of neutron and a-particle emis- 
sions are considered by 

Psv • — (f Pn.evap)^^ Pa,evap\ 

For the application of the evaporation prescription, see our preceding research sum- 
marized in Ref.'^ 

2.3. Treatment of the Fissibility 

First, a long-lived compound nucleus with a certain excitation energy is formed. 
The fissibility (fission probability) is obtained by additional TDDFT calculations 
(cf. the method shown in Ref.l^, where no initial velocity is given. Here we are 
interested in the possibility of fission of a compound nucleus {"^Z): 

^^Zi + ^^Z2, (2) 

where Ai and Zi , which satisfy A = Ai+ A2 and Z ~ Zi+ Z2, are the mass num- 
bers and the proton numbers of the fragments, respectively. First, a configuration 
of the two nuclei {'^^ Zi and "^^^2) at a distance i?o is defined as an initial state 
for additional calculations. Second, choose i?o systematically such that compound 
nuclei with several excitation energies are examined. Note here that the TDDFT is 
a theory in which the total energy is strictly conserved, so that the total energy is 
conserved during the additional TDDFT calculations. Third, the initial many-body 
wave function, which is given as a single Slater determinant, consists of single wave 
functions of two different initial nuclei, where a set of single wave functions are 
orthogonalized before starting TDDFT calculations (cf. the Gram-Schmidt orthog- 
onalization method). 
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Fig. 1. (Colour online) Fissility of compound nuclei based on the liquid drop model (Eq. JSjl). 
Fissilities of ground states are calculated for 0-neutron, 2-neutron, 4-neutron emission channels of 
^**Cm -I- **Ca and for the daughter isotopes of the respective decay chains. 



Fusion, fission and elastic scattering (never touched throughout the time evolu- 
tion) can take place depending on the initial distance (i?o)j because the excitation 
energy is associated with the initial distance. In this way we obtain the possibility 
of fission as a function of the initial distance, and therefore as a function of the 
initial excitation energy. This method brings about the microscopic derivation of 
the fissibility. Note that the present method suggests the fissibility as purely due to 
the collective motion in which the nuclear and Coulomb forces are fully taken into 
account as well as the shell effects and the deformation. 



3. Fissibility of compound nuclei 

3.1. Fissibility based on the liquid drop model 

For comparison, let us introduce the phcnomcnological fissilitj^ (0 < x < 1): 

X-' ,3) 

1; Z^/A > 50, 

where the critical value {Z'^/A)c is roughly equal to 50. In this equation the bal- 
ance between the surface and Coulomb effects is taken into account within the 
liquid drop model description. The value x shows the stability against quadrupole 
deformation (at least); a nucleus can split up into two fragments if x is close to 
1. Note that the fissility is a technical term generally used for the value of x, but 
it exactly means the fissibility. Since this value expresses the instability against 
the quadrupole deformation, it is mostly related with the possibility of symmetric 
fission. For N/Z = 1 nuclei, A — 200 is the threshold value whether x < 1 or not. 



December 4, 2012 2:40 WSPC/INSTRUCTION FILE iwataexon-submit 



Fissibility of compound nuclei 5 



Table 1. Values related with the possibility of fission. Symmetric fission events are considered 
in the former two cases. Asymmetric fission events related with the compound nucleus ^^^Lv 
(no neutron emission channel) are considered in the latter three cases. The contact distance is 
calculated in the liquid drop model (rgA^^^) taking into account the initial deformation of the 
TDDFT calculations. The surface area of the resulting binary system, which is normalized by 
the value for ^''^^Ce + ^*^Ce, is also calculated in the liquid drop model. 



Fission reaction 




Zl X Z2 


Contact distance [fm] 


Surface area (normalized) 


296 Lv ^ 148 Ce 


+ Ce 


3364 


16.3 


1 


284 ^ 142^3 


+ 142Cs 


3025 


16.1 


0.97 


296Lv ^ «Ca- 




1920 


13.7 


0.93 


296l^ ^ 16o + 


280Hs 


864 


12.7 


0.88 


296lv ^ 4He + 


292 pj 


228 


11.5 


0.83 



X = 116^/(296 X 50) = 0.909 in case of the Compound nucleus ^^^Lv. The fissihties 
of ground states for ^^^Lv, ^^"'Lv, ^^^Lv and then several daughter nuclei (from 
different alpha-decay chains) are shown in Fig. 1 in which we found monotonous 
trends depending on both, neutron and alpha emissions. In this paper, instead of 
the fissility x, we present a method of calculating the fissibility in a microscopic 
way. 

3.2. Fissibility based on TDDFT calculations 

Before moving on to the calculation results, the two values related with the fissibility 
are shown for some cases treated in this paper (Table 1). The value Zi x Z2 is related 
with the instability due to the Coulomb energy, while the contact distance is related 
with the threshold energy, depending on whether the nuclear force is operational 
or not. In the case of Eq. ([3|), only the first value is taken into account by assuming 
the operation of the nuclear force. 

Time-dependent density functional calculations with a Skyrme interaction 
(SLy(P|) are carried out in a spatial box of 48 x 48 x 24 fm'^ with periodic boundary 
condition. The unit spatial spacing and the unit time spacing are fixed to 1.0 fm 
and 2/3 X 10~^^s, respectively. The initial states, which correspond to the ground 
states, are prepared by static density functional calculations. 

In order to show a procedure of calculating the fissibility, we have several ex- 
amples. At first, the possibility of symmetric fission is investigated for compound 
nuclei included in an alpha-decay chain. Here we consider the 0-neutron emission 
channel: 

248Cm + 4«Ca ^ 296lv 

^ ^^^F\ + a ^^^Cn + 2a ^84^^ _^ 3^ ^ ... 

as a limiting case. More concretely let us study the symmetric fissibilities of 
^^^Lv and ^^^Ds corresponding to the cases with -^^Zi = "^^22 = ^^^Ce and 
■^^Zi = ^^^^2 — ^^^Cs (in Eq. ([2])), respectively. Next, the possibility of asym- 
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Fig. 2. (Colour online) Systematic TDDFT results for two compound nuclei. Results starting with 
different binary systems are shown for different initial distances _Ro = 10, 12, 14, 16, 18 fm. The 
initial compound nucleus is set to ^^^Lv except for a symmetric fission: ^**Ds — > ^''^Cs + ^'*^Cs. 



metric fission is investigated for tire compound nucleus ^^^Lv. Although the fissility 
X is the value for symmetric fissions, the fissibility can be considered using the 
present method. The following fission events arc considered: 



296Lv ^ 


Cm 


+ 4"Ca 


296Lv ^ 


280jjg ^ 


- 1^0, 


296Lv -> 


292pi _j_ 


^He, 



where the mass asymmetries are 0.16, 0.06 and 0.01, respectively. 

The systematic TDDFT results with different Rq values are summarized in 
Fig. 2. It is notable that no fission appears in all asymmetric cases, where the 
attraction due to the nuclear force is considerably larger than the repulsion due to 
the Coulomb force. It means that the emission of small nuclei such as alpha-particles 
(alpha-decay) cannot be explained by collective dynamics in which the competition 
between the surface tension and the Coulomb force are fully taken into account. 
In this context, another effect such as the thermal instability and the resulting 
cooling effeclPl might play a role in the emission of small nuclei, instead. It is also 
suggested that there is a transition of fission mechanism of compound nuclei; from 
the mechanism mostly due to collective dynamics to another mechanism. According 
to Table 1, the Coulomb energy of symmetric fission (min(^i/^2, •^2/^1) = 1) 
is 1.75 times larger than that of asymmetric fission with min(^i/^2, A2/A1) = 
0.16, while the surface energy ratio between them is only 1.08. It implies that the 
transition of fission mechanism arises from the different mass dependence between 
the surface and Coulomb effects. This transition is identified by the threshold value 
for the mass asymmetry ratio. 

For the results of symmetric fissions shown in Fig. 2, no significant difference 
is noticed between the two cases. In those calculations the incrementation of Rq is 
fixed to 2 fm, which is too large to identify the difference between the two cases. On 
the other hand, difference in the results of asymmetric fissions shown in Fig. 2 can 
be essentially understood by the difference of the contact distance (see Table 1). 

For each reaction, let us define the fissibility (effective fission probability) as 
follows: 



(4) 
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Fig. 3. (Colour online) Fissibility of compound nuclei. 



: for very mass-asymmetric cases, 
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intermediate mass-asymmetry region; 0.16 < m'm{Ai/A2, A2/A1) < 1. 



where ajig is the normahzcd cross-section satisfying X]fio<J?c ~ ^ ^^'-^ ^^'^ 
contact distance), and Pfis,Ro = 1,0 is the probabihty of fission; Pfis,Ra — 1 when 
the fission is calculated for the given initial distance i?o- Therefore is the cross- 
section for realizing the initial state of TDDFT calculations with the initial distance 
i?o. This parameter (0 < O < 1), which is defined for each reaction (Eq. ^) with 
fixed -^Z, -^^ Zi and ^^^2, is the microscopic correspondence of the fissility x. O 
is equal to if fission never appears, and to 1 if fission necessarily appears. It is 
notable that, as is seen in the symmetric cases shown in Fig. 2, the fissibility can be 
sometimes equal to zero using this fissibility (cf. fissility x). In order to clarify the 
difference depending on the mass asymmetry, let us introduce a simplified fissibility 
following Eq. (|4]). 

Csimp = max(P/js,_Ro). (5) 

This parameter (Osimp = 0, 1) identifies if fission possibly appears (Osimp = 1) 
due to the collective dynamics (essentially the competition between the surface 
and Coulomb effects) or not (Osimp = 0). Figure 3 shows the behaviour of the 
simplified fissibility within the performed calculations, where additional calculations 
with Rq = 0, 2, 4, 6, 8 fm are made for three mass asymmetric cases in order 
to determine the value of Osimp- Such a mass dependence cannot be described 
by the fissility x. The investigation of the intermediate mass-asymmetry region 
(0.16 < min(^i/^2, A2/A1) < 1) will be studied in the future. 

4. Summary 

Compound nucleus fission is a complicated scenario in which the collective and sin- 
gle particle degrees of freedom, thermal effects (e.g., cooling effect leading to the 
localization), the collectivity due to the pairing, and other effects are activated on 
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a casc-by-casc basis. Among these factors, much attention has been paid to the 
coUectivity in this paper. A procedure of calculating the fissibility of compound 
nuclei has been presented. The present method allows us to have a microscopic 
treatment for the fissibility for both symmetric and asymmetric events. In partic- 
ular, using this method, the shell effect, the deformation of the fission fragments 
and the fermionic property (anti-symmetrization of many-body wave function) are 
introduced without having any additional treatments. 

For testing, the method has been applied to analyse the fissibility of the com- 
pound nucleus ^^^Lv. The difference between symmetric and very asymmetric fis- 
sions has been clarified. In fact it turned out that very mass-asymmetric fissions 
cannot be explained by the collective dynamics, although the competition between 
surface and Coulomb effects can be fully taken into account. 
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